Higher-order Kerr terms allow ionization-free filamentation in gases 
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Abstract 

We show that higher-order nonlinear indices (77-4, ri6, n\o) provide the main defocusing con- 
tribution to self-channeling of ultrashort laser pulses in air and Argon at 800 nm, in contrast with 
the previously accepted mechanism of filamentation where plasma was considered as the dominant 
defocusing process. Their consideration allows to reproduce experimentally observed intensities and 
plasma densities in self- guided filaments. 
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The filament at ion of ultrashort laser pulses in gases [1] attracted a lot of interest in the 
last years because of its physical interest as well as its potential applications (2- 5] . Filaments 
are self- channeled structures propagating over many Rayleigh lengths without diffraction. 
They are generally considered to stem from a dynamic balance between Kerr focusing and 
defocusing by the plasma generated at the non-linear focus. Numerical simulations based on 
this balance report a core intensity of several 10 13 W/cm 2 and typical electron densities of 
several 10 16 cm -3 3I, U]. Consequently, plasma ionization is generally admitted as necessary 
for an ultrashort pulse to experience self- channeling in gases. 

But the plasma density provided by this description of filamentation appears overesti- 
mated as compared with experimental measurements. As reviewed in [6], such measurements 
are dispersed over several orders of magnitude, especially due to different focusing conditions 
and divergent assumptions about the core diameter of the filaments, but the electron density 
in a filament generated by a slightly focused beam is more likely to amount to 10 14 — 10 15 
cm -3 6]. This value, as well as the discrepancy by more than one order of magnitude with 
numerical simulations, was recently confirmed [7||. The observation of so-called plasma- free 
filamentation (si, lofl, as well as the consideration that a balance between the instantaneous 
Kerr term and the time-integrated plasma contribution implies strongly asymmetric pulse 
shapes [10|, periodically led to challenge the role of plasma in laser filamentation. 

However, up to now, no other process seriously challenged plasma as the main defocusing 
process balancing the Kerr self-focusing. Nurhuda et al. proposed that the saturation of 
the nonlinear susceptibility x^ should be taken into account Such saturation can be 
described as negative higher-order Kerr terms. The nonlinear index of air induced by high- 
power femtosecond laser pulses can be written as An Ker r = n 2 /+n4/ 2 +n 6 / 3 +n 8 / 4 +... , where 
/ is the incident intensity and the n 2 *j coefficients are related to y( 2 *j +1 ) susceptibilities. This 
nonlinear index is generally truncated after its first term, n 2 |2rl5(, mostly because of the 
lack of data about the values of the subsequent terms. 

Numerical works have investigated the influence of the quintic nonlinear response on 



the propagation dynamics in gases, although without knowledge of its value |12Ml6j. They 
showed that 714 is negative, i.e. the x^ susceptibility is a defocusing term. It tends to 
stabilize the propagation of ultrashort laser pulses in air and to decrease both the electron 
density and the maximal on-axis intensity. Consequently, the losses due to multiphoton 
absorption (MPA), which lead to the end of the filamentation, are reduced and pulse self- 
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channeling is sustained over longer distances. However, plasma generation still appeared 
as necessary for filament stabilization. Moreover, the value of 77,4 was set arbitrarily, which 
limits the conclusiveness of these studies. Finally, the lack of data prevented any evaluation 
of a possible effect of the further-order nonlinear refractive indices. 

However, the higher-order Kerr indices have recently been measured in N 2 , O2 and Ar 
by Loriot et al. [17]. The reader is referred to this work for a detailed description of 
this experimental determination. In this Letter, we investigate their influence on numerical 
simulations of laser filamentation. We show that their values are sufficient to provide the 
dominant contribution to the defocusing terms of self-channeling. Their implementation in 
numerical simulations yields the experimentally observed plasma density. As a consequence, 
contrary to previously held beliefs, a plasma is not required for the observation of filamen- 
tation. Rather, plasma generation can be considered as a by-product of the self-guiding of 
laser filaments. 

We implemented these nonlinear coefficients into a numerical model describing the propa- 
gation of ultrashort high power pulses . We consider a linearly polarized incident electric 
field at Aq=800 nm with cylindrical symmetry around the propagation axis z. The scalar 
envelope e(r, £, z\ assumed to vary slowly in time and along z, evolves according to the 
propagation equation: 



where /co=27m /Ao and uj = 2ttc/Xq are the wavenumber and the angular frequency of the 
carrier wave respectively, n is the linear refractive index at A , k" = f^l^o is the second or- 
der dispersion coefficient, p at the neutral atoms density, p the electron density, p c = e^mujl/e 2 
is the critical electron density, m being the electron mass and e its charge. Vl^(|6:| 2 ) and 07 
are the photoionization probability and the inverse Bremsstrahlung cross-section of species I 
respectively (with ionization potential £/j), and t refers to the retarded time in the reference 
frame of the pulse. The right-hand terms of Eq.(pQ) account for spatial diffraction, second 
order group- velocity dispersion (GVD), instantaneous nonlinear effects (i.e. the nonlinear 
refractive index of air, up to the n 8 term), plasma defocusing, inverse Bremsstrahlung and 
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Species 


cm 2 /W) 


cm 4 /W 2 ) 


cm 6 /W 3 ) 


cm s /W 4 ) 


N 2 


1.1±0.2 


-0.5±0.27 


1.4±0.15 


-0.44±0.04 


o 2 


1.6±0.35 


-5.2±0.5 


4.8±0.5 


-2.1±0.14 


Air 


1.2±0.23 


-1.5±0.3 


2.1±0.2 


-0.8±0.06 



Table I: Coefficients of the nonlinear refractive index expansion of N2 and O2 at 1 bar pressure, 
and interpolation to air, as used in the present work 

multiphoton absorption respectively. As compared with previously published data [17], we 
used values of the higher-order refractive indices (Table [I]) incorporating the correction for 
the coherent artifact [191 ]. i.e. adequately substracting its electronic contribution at play in 
the original measurement of Ref . (l7| • This correction results in dividing each n 2 *j term by 
j + l. Owing to the short pulse duration (30 fs) used in the simulations, the delayed orienta- 
tional response is disregarded. The propagation dynamics of the electric field is coupled with 
the density of the electrons originating from the ionization of both O2 and N 2 : p = Po 2 +Pn 2 - 
This density is governed by the muti-species generalized Keldysh-PPT (Perelomov, Popov, 
Terent'ev) formulation [3, (|. 

We used this model to simulate the propagation of an ultrashort pulse typical of 
laboratory-scale experiments: 1 mJ energy, 30 fs FWHM pulse duration without initial 
chirp (hence, about 3.9 critical powers Per), an initial waist of = 4 mm, a focal length 
/ = 1 m and a pressure of 1 bar. Figures [1] and [2] compare the numerical results of the full 
model implementing Kerr terms up to n% and of the classical model, where the Kerr term 
is truncated to n 2 . Both models lead to self-guided filaments. The full model yields a lower 
maximum intensity (31.6 TW/cm 2 vs. 78 TW/cm 2 ), although these values lie within the 
range of published experimental data in comparable conditions ^2-5\. On the other hand, 
the full model predicts an electron density 40 times below the classical one (1.1 • 10 15 cm -3 
vs. 4.2 • 10 16 cm -3 ). While the latter value is comparable with the output of other numeri- 
cal works the full model agrees with the available experimental measurements of the 
electron density jf], 7]. 

Note that, with the considered parameters, the full model yields a more strict intensity 
clamping than the classical one J20I]. It predicts an intensity constant within 20% over 15 cm 
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— Full model — Classical model ■■■ Full model without plasma 




Figure 1: (a) On- axis intensity and (b) Plasma density as a function of the propagation distance 
for the classical model (considering only ri2 term of the Kerr index and the plasma defocusing), the 
full model, as well as the full model without plasma. 
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Figure 2: Fluence distribution in air as a function of the propagation distance for the full model 
(a) and the classical model including 712, ionization and GVD only (b). The white lines display the 
quadratic radius as a function of the propagation distance. 



vs 9.5 cm in the case of the classical model), a length well comparable to experimental data 



reported to date in air for mJ-pulses 
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221]. This stricter clamping can be explained 



by the lower electron density, which results in weaker multiphoton losses, allowing a slower 
decay of the filament intensity and ionization. The full model also yields a narrower output 
spectrum (Figure [3]), which better fits experimental data in air ^, U]. It should therefore be 
considered as the reference model for numerical simulations of filamentation. Note that the 
almost symmetric shape of the spectrum is due to the neglection of self-steepening 
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Figure 3: Spectrum after 2 m propagation in air at atmospheric pressure. 

On the other hand, neglecting the ionization in the full model (see Fig. [D(a)) almost does 
not affect the simulation output. This shows that, in contrast to the classical understanding 
of filamentation in gases, the self-guiding process and plasma generation are almost decou- 
pled. Instead, the negative higher-order nonlinear indices 714 and n% constitute the dominant 
regularization terms leading to filamentation in air at atmospheric pressure. This limited in- 
fluence of the ionization on the filamentation dynamics when higher-order non-linear indices 
are adequately considered sheds a new light on the possibility of ionization-free filamenta- 
tion [8], which appears as a natural possibility in the context of the full model. Still, the 
dominant contribution of higher-order Kerr terms does not prevent ionization (Fig. [lb), 
which may contribute e.g. to the conical emission. 

We checked that the above conclusions are not restricted to a particular set of values of the 
non-linear refractive indices. Indeed, qualitatively comparable results have been obtained 
when varying the indices by several tens of percent, comparable with the experimental 
uncertainties on the non-linear indices. Furthermore, to compare the above molecular results 
with an atomic gas where no molecular orientation occurs, we performed simulations for 
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Table II: Coefficients of the nonlinear refractive index expansion of Ar at 1 bar pressure, as used in 



the present work 
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Figure 4: (a) On-axis intensity and (b) Plasma density as a function of the propagation distance 
for the classical model (considering only ri2 term of the Kerr index and the plasma defocusing) and 
the full model, in Argon under 1 bar pressure. 

Argon, where the ionization potential is close to that of the air molecules [23], thus behaving 
in a similar manner as far as ionization is concerned. As in the case of air, we refined the 
corresponding indices to take the coherent artifact into account. The resulting values are 
summarized in Table [III Like in air, the full model yields lower filament intensity (28.5 
TW/cm 2 vs. 80.9 TW/cm 2 ) and electron density (5.2 • 10 13 cm -3 vs. 4.1 • 10 16 cm -3 ) than 
the classical model (Figured]). Also, the evolution of the fluence profile as a function of 
propagation distance (Figure [5]) is quite similar in both models. 

The space-time dynamics shows more differences between the full and the classical models 
(Figure Ej). In both cases, the pulse splits into two sub-pulses around 1.05 m propagation, 
but the full model predicts an almost symmetrical temporal profiles pattern all along prop- 
agation, while the classical model yields a largely asymmetric one. This behavior illustrates 
the different temporal dynamics of higher-order Kerr terms as compared with the plasma 
generation. The former is an instantaneous phenomenon depending only on the intensity. In 
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Figure 5: Fluence distribution in 1 bar Argon as a function of the propagation distance for the full 
model (a), and the classical model including n2, ionization and GVD only (b). 

contrast, the plasma generated during the pulse accumulates, resulting in an ever growing 
contribution. As a consequence, the leading edge of the pulse propagates in a low plasma 
density while the trailing edge is more defocused by the much higher electron concentration 
it encounters. Moreover, the lower losses due to the lower plasma density in the full model 
allows a slight refocusing cycle around 1.15 m, which is not predicted by the classical model. 
The results of the full model stay unaffected when the plasma is not taken into account (e.g 
the peak intensity only increases by 0.6 %), which confirms that the filamentation process, 
including the pulse splitting is indeed driven by the higher order Kerr terms when they are 
considered. 

These differences in time-space dynamics illustrate the interest of implementing all orders 
of the Kerr effect in numerical simulations of filamentation in gases. Since successive terms 
ri2*jP of the Kerr index are of alternate signs and have comparable values at an intensity of 
about 30 — 35 TW/cm 2 [Tjj, the inclusion of all terms up to n 8 in air (resp. in Argon) 
is necessary to adequately simulate the propagation of filamenting ultrashort pulses. 

The observation that ionization, as well as GVD, almost do not affect the results of the 
full model provides an opportunity to speed up the numerical simulations. Neglecting the 
ionization typically cuts the computation time by a factor of 3 with little impact on the 
result in the conditions shown above. A parametric study would be necessary to determine 
the conditions, and especially the wavelengths and materials where such approximation is 
legitimate. Such study shall compare the intensities yielding a dynamic balance of the Kerr 
terms on one side, and between Kerr and plasma contributions on the other side. In air, 
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Figure 6: Space-time dynamics of filament ation in 1 bar Argon for the full model (a), and the 
classical model including 77,2, ionization and GVD only (b). Both models yield pulse splitting around 
1.05 m propagation distance, but the full model where fllamentation is driven by the instantaneous 
Kerr effect results in a more symmetrical temporal dynamics. 

where these intensities amount to 31.6 TW/cm 2 and ^ 78 TW/cm 2 , respectively, the lower 
intensity for pure Kerr balance ensures the domination of the latter process. Depending 
on the respective values of the higher-order non-linear indexes and ionization rates, the 
respective balance intensities may switch, leading to the domination of the Kerr-plasma 
balance. 

In conclusion, we have shown that the recently measured higher-order nonlinear indices of 
air (up to n$) or Argon (up to 77,10) dominate both the focusing and defocusing terms implied 
in the self-guiding of ultrashort laser pulses in these gases. As a consequence, contrary to 
previously held beliefs, a plasma is not required to generate fllamentation in gases, and its 
generation is quite decoupled from the self-guiding process. Instead, fllamentation is, at least 
in the considered conditions, governed by higher-order nonlinear indices. The usual definition 
of a filament as a dynamic balance between the n 2 Kerr self-focusing and defocusing on the 
plasma shall therefore be revisited. Filament ation in gases rather appears as a nonlinear self- 
guided propagation regime sustained by a dynamic balance between nonlinear self-focusing 
and defocusing effects. Depending on experimental conditions, the latter can include higher- 
order Kerr terms and free electron with respective weights depending on the propagation 
medium considered. 
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